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ABSTRACT 


^  Prcscincad  in  this  report  ar«  toLdl  aud  jii>i:rur'>^-itlon  hj.; 
degree  of  freedom  equations  of  muclon  and  tneir  c  rlvations. 

Magnus  force  and  moment  were  considered  ntii;ii>^.i oJe  and  not 
included  in  this  derivation,  however,  when  desired  they  can 
be  incorporated  in  the  existing  models.  Definitions  oi  control 
force,  gravity,  thrust  and  damping  tarms  sire  presented  in  general 
form  so  that  various  control  schemas  and  particular  functions  of 
tha  terms  can  ba  utilised  without  estoneive  modification  to  the 
equatlone  of  motion.  Theee  equations  may  be  prograimud  on  a  large 
scale  digital  or  analog  computer. 

The  majority  of  tha  work  presented  herein  wos  developed  for 
previous  programs  and  adapted  for  use  in  the  Penetration  Aids  Program 
under  Air  Force  Contract  AF04  io VA ; >25 . 


OBJECT; 


The  object  or  this  report  i;.  Lo  aorivo  und  preuoni,  .  ix-tio.p'tfe-Oi 
equetlODS  of  notion  for  •  ri^id  ba.lliritic 

nfmODUCTION; 

Fli^t  nechanics  problems  involving  non-aplnnln^  vehicles  c«n  be  solved, 
in  genenlf  vlth  dynanloel  equations  which  are  independent  of  the  vehicle's 
point  aass  trajectory,  in  the  spinning  body  problem,  however,  the  orientation 
of  the  vehicle  and  its  motion  are  mutually  related  as  a  result  of  aurodynamic 
and  inertial  coupling  and,  therefore,  must  be  studied  with  a  six-degree -oi'-freedoa 
model.  This  report  presents  the  derivation  of  several  s Lx-degree-of-freedom 
models  for  a  general  problem.  It  has  been  assumed  in  this  derivation  that  the 
earth  la  spherical  and  non-rotating  und  that  body  spin  rater,  are  small  enough 
to  neglect  Magnus  forces. 

nie  c.g.  fixed  part  of  the  motion  of  the  body  is  expressed  non-holonomically, 
In  equations  (3)i  as  dynamical  equilibrium  aoout  any  three  riglit -handed, 
orthogonal  cartesian  axes  fixed  in  the  body.  Tlte  coordinates,  p,q,  and  r,  of 
equations  (3)  ere  known  classically  uv  the  rectilinear,  but  ncixa^oionuoilc 
coordinates  of  rigid  body  rotation  dynamics.  These  coordinates  can  be  expressed 
in  terms  of  holonooiic,  but  curvilinear  coordinates  such  as  Euleriun  Angles. 
Equations  (U)  represent  transformation  oi  p,q,  and  r  to  tlie  particular  Kulerian 
Angles  f  ^  t  end  By  means  oi  truns  format  ion  (4),  the  dynamical 

equilibrium  expressed  non-holonomically  by  equations  (3)  can  be  rirexpressed 
holonoiaically  as  it  is  in  equations  (!?).  The  dynamical  e(ii!Qtiont  lor  body  rotation 
are  non-holonomic  In  p,qi  and  r  because  tne  inertia  terms  in  (3)  are  only  once 
Integrable  without  the  adjunction  of  equations  (4). 


i 


Immodlately  subsequent  to  equations  (4),  in  this  report,  the  question  of 
what  coordinates  to  use  io  decided.  Two  approaches,  designated  moasnt  derivation 


A  and  nooent  derivation  B,  are  adopted.  The  starting  points  for  these  are 
equations  (S)  for  A  and  equations  (3)  for  B.  In  both  "derivations",  the  applied 
torque  terms  are  expressed  in  terms  of  Bulerian  Angles,  ^  ^  ,  and  . 

Only  in  A,  however,  are  the  inertia  terms  also  expressed  in  if,  t  0  ,  and  ^  . 

An  Interim  result  in  A  is  equations  (18)  which  are  a  direct  combination  of  the 
inertia  terms  of  equations  (^)  with  the  applied  torque  terms  of  equations  (17). 
Equations  (19)  through  (23)  are  concerned  with  the  translational  equations 
of  motion  of  the  body  in  the  space-fixed  coordinates,  x,  y,  and  z  of  translation 
and  in  the  space-fixed  coordinates  7^  ,  0  ,  und  ^  of  orientation.  Equations 
(16)  together  with  equations  (23)  constitute  a  complete  set  of  holonomic,  but 
non-linear,  dynamic,  six-degree -of-freedom  equations  of  motion  of  the  rigid  body 
in  the  all -space -fixed  coordinates  x,  y,  z,  ^  ,  &  ,  and  ^  .  Because  the 
individual  rotation  equations  of  this  system  represent  dynamical  equilibrium  of 
cosqpoixents  along  body-fixed  axes,  only  constant -valued  second  and  product  moments 
of  inertia  occur.  Equations  (18)  plus  equatlous  (23),  as  a  system,  are  incomplete 
in  that  flexible -body  degrees -of-freedom  are  excluded  and  are  inexact  In  that 
(1)  magnus  forces  81*0  neglected  and  (2)  a  spherical,  non-rotating  earth  is  assumed. 

Further  development  under  mi  ment  derivation  A  is  desirable  because  bodies 
either  symmetric  about  the  roll  axis  or  nearly  symmetric  about  the  roll  axis  are 
a  frequent  object  of  dynamical  analysis.  In  the  translational  degrees -of-freedom, 
bodies  with  roll  symmetry  can  be  governed  no  more  simply  than  by  equations  (such 
as  equations  (25)  which  govern  unsymaetrical  bodies.  In  the  rotational  degrees- 
of-freedom,  however,  a  form  simpler  than  equations  (I6)  is  possible  for  the  roll 
symmetric  body.  Such  symplificatlcn  obtains  by  the  following  reasoning.  Instead 
of  equations  (I8)  which  represe^nt  d/namic  equilibrium  of  components  of  applied 
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torqut  along  body-flxad  axas  vlth  cooqponantB  of  t;->«  applioabl*|  rtaetlva,  inortla 
torquaa  along  tba  oaan  axas,  aquatlona  i^ill  based  on  ^  ,  <9  |  and  but 
•pacifying  dynanical  aqulllbrlun  of  corresponding  conponants  along  any  suit¬ 
ably  United  aat  of  three  axes  can  be  used  with  equal  efficacy  to  fix  the 
rotational  behavior  of  the  unsysnetrical  body.  In  general,  a  (pre -solution) 
dioloe  of  conponent  axes  flxed-ln-space  tends  to  produce  an  inertiaily  uncoupled 
set  of  equations,  but  also  leads  to  variable  second  and  product  nonsnt  of  Inertia 
coeffleiente.  Conatant-lnertla-coeffielent  equations  applicable  to  the  unsynsstrlc 
body  require  a  choice  of  conponent  axes,  such  as  that  choice  resulting  In  equations 
(18),  iriilch  follow  the  body  In  yav,  pitch,  and  roll.  For  bodies  synsetrlcal 
to  the  roll  axis,  a  choice  of  conponent  axes  which  follow  the  body  only  In  yav 
and  plt^  gives  constant-inertla-coefflclent  equations  with  a  fora  less  inertiaily 
coupled  than  equations  (IC*;.  Obe  x',  y',  and  z*  axes  system  follows  the  body  in 
yaw,  pitdi,  and  roll  diile  the  x^^,  y^^,  and  z^^  axes  system  follows  the  body  only 
in  yav  and  pitch.  Bquations  (S6)  express,  for  any  static  moment  vector,  components 
along  the  Instantaneous  x^^,  y^^,  and  z^^  axes  in  terms  of  components  along  the 
Instantaneous  x',  y',  and  z*  axes.  Applicatlcxi  of  equations  (26)  to  the  inertia 
torque  terms  and  to  the  applied  torque  terms  of  equations  (10)  yields  equation 
of  dynamic  equilibrium  In  rotation  about  the  ,  y^^  ,  and  z^^  axes.  These 
equations  are  valid  for  the  unsymmetric  body.  Subsequent  specialization  to  the 
roll  sysnetric  body  is  readily  accomplished  by  setting  every  second  moiBent  of 
Inertia  about  any  axis  perpendicular  to  the  roll  axis  equal  to  a  single  value,  I, 
and  by  setting  all  product  moments  of  Inertia  equal  to  zero.  The  result  is 
equations  (27)  which  are  exact  for  the  roll  syimistrlc  body  in  the  same  sense  that 
equations  (16)  are  exact  for  the  unsymmetric  body.  Equations  (2^)  together  with 
equations  (2?)  constitute  a  coaplete  set  of  holonomlc,  but  non-linear,  dynamic, 
six -degree -of -freedom  equations  of  motion  of  the  roU-symnetrlc,  rigid  body  in 
the  all-space-fixed  coordinates  x,  y,  z,  ^  ,  9  ,  and  .  Because  the 

individual  rotation  equations  of  this  system  represent  dynamical  equilibrium  of 
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conrponents  along  axes  which  follow  th«  body  in  yaw  and  pitch,  and  because  they 
are  based  on  the  assumption  of  roll -symmetry  of  the  body,  only  constant  valued 
second  moments  of  inertia  occur  and  product  moments  of  inertia  do  not  occur  At  all. 

Continued  development  under  moment  derivation  A  consists  of  forming  pertur> 
bation  equations  of  motion  from  the  system  composed  of  equations  (25)  together 
with  equations  (27).  This  is  accomplished  by  subtracting  equations  corresponding 
to  a  standard,  zero  angle  of  attack  trajectory  from  equations  (25)  and  (27), 

FOr  this  standard  trajectory,  equations  (25)  and  (2?)  reduce  to  equations  (28), 
When  equations  (28)  are  subtracted,  then,  from  equations  (25)  and  (27),  the 
result  is  equations  (29)  which  are  the  desired  perturbation  equations  of  motion. 

For  applications  to  roll-symmetric  bodies  and  where  the  yaw  angle,  ^  ,  and 
the  miscellaneous  angles  and  small  enough  to  linearize, 

equations  (29)  reduce  to  the  (approximate)  equations  (30).  Equations  (31) 
through  (34)  are  concerned  with  introducing  the  parameters  and  Oty  and 

with  further  linearization  of  equations  (30),  Equations  (34)  are  the  final 
result  under  moment  derivation  A  and  represent  approximate  equations  of  motion 
for  the  roll  symmetric  body.  Assumptions  incorporated  in  equations  (34)  and 

•  •  s  * 

A  - 


not  Incorporated  in  equations  (30)  are  that  (  ^ ), 


and 


V 


are  small  enough  to  linearize. 


under  moment  derivation  B,  only  the  rotation  equations  are  treated.  Hie 
translational  equations  are  the  same  for  B  as  for  A,  Equations  (35)  are  a  set 
of  three,  non-holonomlc,  six -degree -of -freedom  equations  for  dynamic  equilibrium 
of  torque  components  along  the  x',  y',  and  z'  axes  system  in  nine  coordinates, 
p,  q,  r,  X,  y,  z,  ,  and  ,  Together  with  the  three  translational 

equations  (25),  and  with  the  three  "constraint”  equations  (4),  equations  (35) 
form  a  complete  set  of  nine,  non-holonomic,  six-degree-of-freedom,  equations  of 


4. 
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■otion  Applicable  to  th«  gtooralf  unsyBBotriCf  rigid  body.  Equation  (3?)  are 
foraed  by  setting  the  applied  torque  terns  of  equations  (17)  equal  to  the  reactive, 
inartia  torque  terns  of  equations  (3)  and  represent  dynanlc,  rotational  equilibrium 
about  axes  vhleh  follow  the  body  in  yaw,  pitch,  and  roll. 

For  the  roll  synnetrlc  body  in  notion  linearizable  with  respect  to  &  ,  /^p  , 
A.  ^  ^  ,  and  ,  equations  (36)  supplant 

equations  (35)  for  expressing  rotational  equilibrium  and  equations  (3^)>  (3^b), 
and  (3^c)  supplant  equations  (2^)  for  translational  equilibrium.  Equations  (36)  differ 
from  34d,  3Ve,  and  3^f,  not  only  because  of  the  use  of  p,  q,  and  r  in  expressing 
the  reactive  inertia  terns,  but  also  because  equations  (36)  represent  dynamic 
equilibrium  about  axes  which  follow  the  body  in  yaw,  pitch,  and  roll  ^lle 
equations  3bd,  3^e,  and  3^f  represent  dynamic  equilibrium  about  axes  ^Ich  follow 
the  body  only  in  pitch  and  yaw. 


COHCLUBIdlS: 

It  is  concluded  that  this  report  presents  useful  forms  for  the  following 
sets  of  equatlcxus. 

Eqns.  (16)  Total  rotational  equations  in  terms  of  Eulerlan  angles. 

Eqns.  (25)  Total  translational  equations  along  space  fixed  axes. 

Eqps.  (29)  Sisqplified  six -degree -of -freedom  perturbation  equations. 
Sqna.  (34)  Linearised  six-degree-cf^freedom  perturbation  equations, 
tiqns.  (35)  Total  body  moomnt  equatlcxis  in  terms  of  body  angular  rates. 
E^is.  (36)  Linearised  bod:  moment  perturbation  equations  in  terms 
of  body  angular  rates. 


Derivation  of  the  Equations  of  Motion 

The  translational  equations  of  motion  ore  written  with  respect  to  a 
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spaced  fixed  axes  system  In  order  to  measure  true  accelerations  In  fixed 
space  to  which  Newtonian  lawj  apply.  The  moment  equations  are  written 
relative  to  body  fixed  axes  with  origin  at  the  center  of  gravity  and 
parallel  to  the  axes  of  symmetry  of  the  body.  Body  fixed  axes  eliminate 
the. need  for  considering  the  moments  of  Ineftle  as  a  function  of  vehicle 
orientation.  The  thrust,  control  force,  gravity,  and  daiiqping  definitions 
ara  general  since  they  will  be  uniquely  defined  for  each  problem.  The 
moment  or  rotational  equations  of  motion  are  considered  first. 

Section  I  -  Rotational  Equations  of  Motion 

From  reference  1  and  2  the  moments  of  momentum  of  a  missile  about  its 
body  fixed  axes  system  are  defined  as: 

Moments  of  Momentum 


(1) 


Prom  flgiure  1,  the  moments  about  the  body  fixed  axes,  defined  aa  the  rate  of 
change  of  moment  of  momentum  about  those  axes,  are  expressed  as: 


(2) 


^oll  (Body)  .  -  A  h.,  q 

“pitch  (Body)  .  V-  /fy  r* 

“yaw  (Body)  » 


6. 


LfiHjiteWi 
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Substitutinf;;  equations  (1)  Into  equations  (2),  presents  the  body 
moment  equations  in  the  form: 

«Rou  .  i-Cr'-prA'+frCti'i,) 

H».-  .  ri^.  t  Cf-f?/ T,V  + 

In  these  equations  the  terms  cuntuining  the  time  derivatives  of  moment  of 
inertia  have  been  neglected. 

The  angular  rates  about  the  body  axes  in  terms  of  Eulerian  angles  are 
shown  from  figure  2  to  be: 

p :  ^  ^ 

(*♦)  q.  S 

^  -  9  ^4^^^ 

At  this  point  in  the  derivation  a  choice  must  be  made  to  proceed  by 
either  expressing  body  moment  In  terms  of  Eulerian  angles  and  their  rates 
or  In  terms  of  body  angular  rates  and  accelerations.  Both  are  presented, 
however,  the  former  Is  taken  first  in  this  reports  The  derivations  will  be 
labeled  A  and  B  for  their  reepefctlve  occasion  in  the  report. 

Koment  Derivation  A 

Substituting  the  valuei.)  of  p,  q,  and  r  as  defined  above,  and  their 
derivatives  into  equations  3,  the  body  moment  equations  becomes t 
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Z.  ‘  'Hull  •  4'  ^  ^  ^ 

-A  ^*^^'h'j.t9-^c/  IvLf  A^/'  C*'^ 

Cci&c  t'^ip  0  ■*<,%.  0  ^  —  0  .-*^v  ^  -  4%x.0 

( J'l)  -A  (^U'^  •****•  o.  ^  Cf:i  0  A/l.  ^<r>i  <P  0  cun.  <fi 

CUJL  0 

A  ^  r^0  fcc^!^^-’  —  A©  ^  c^2.  ^  cui.  0  cf> 

7  ^Itch  •  Ce3..<^  »-  y^>y.  ^  -A  ^  i^^V'JL  0  <^^7^0  <MA. 

A^^  ci^-2.6cex.  </>  ) 

+(!,.- I^)(H  e  ¥>-^'  yd*'t\.  Ocux0c:aL^ 

(5b)  A»^^/ ^-A  A  ^  i^Cr^ O C*:t! ^ 

''•^/t^f  jjL u^H.  0  ~  C*^7L  0 ^0 

“A  0 ^f\.^  CUra^J^^ 

Ac^  jk'  a  A  ^’a.  (^  A  A 

—  /^*4X,  0 c^0  WfC'  ^--  ^  C3rx0 Ccn.  (f>) 
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— e  Aam3^  —  e(>  Co-a»^) 

*♦“  (Xy-X^’J^^Co-d.^  ■4-^'^C^B^i^4>’^y^jiiftJ0  C^P  -v 

—  JitV(^b  ‘‘^jduvoe  *  W (hi^ ^ 

-  ^  dyd^P^ 

-  ■^v^‘  (©G>^^  +2(}>V/  Co-u©  eo^4^  -2©4>>4^4> 

-y^  -h  f  Cix9^^ 

c^eAt^P  -  <^*" 

•4-z^y>4^e  “♦" 


Applied  Moment 


The  total  applied  moment  about  a  body  axis  Is  comprised  of  an 
aerodynamic  force  moment,  a  control  force  moment,  a  thrust  vector 
mlaallgnment  moment,  and  a  damping  moment.  This  can  be  written  in  equotion 


form  as: 


(6)  M»  *  Mm,»e  Mr  r  +  M^.  -^ 


Moment  is  expressed  in  vector  analysis  as: 

utting  ~^'A,y' i-^'Am 

and  -»  —  ^ 


where^  is  the  vector  from  the  'DodyC.G.  to  the  application  of  the  aerodynamic  force 
and  where  i*  ,  J* ,  and  k*  are  unit  vectors  alon,^  the  body  y/ ,  y'  and  z* 
ixes, respectively. 

The  above  cross  product  has  the  following  components: 


^ero 


(7)  »y:, 


"  -  ftAnt 

Ay’ 


The  aerodynamic  force  components  as  shown  in  figure  3  a^re: 


*/Uf 

where : 


“‘4 


‘^CAO 


With  the  vector  transformation  on  figure  2,  velocity  can  be  expressed 
along  body  axes  so  that  Fjj  becomes: 

t  .  Qi^-LT'lkC^(bAi^9Co<V-(^'t>^F^  N 

+ Y  (^  ><>  ^  ® ^ 

+  2 

The  y^and  z^coraponents  are: 

*•  4  \v^c4. 

— 2F  C<M1^ 


-.<(Q^(^  ti^e  Ai^f-.^<^ 

-i  Ci^(^Co<y^ 


Kje  monient  arm  components  are  defined  as: 


r^t  m  —  (CG-CP)  vhere  CG  and  CP  are  irveasured  from  an  aft  reference  point 

(9)  ry  s  ^  CPy' 

T^$  «  CPy 


CPy'and^CPj*  include  the  misalignment  of  the  center  of  pressure  and 
the  center  of  gravity. 


Substituting  equations  8  and  9  into  equations  7/  presents  the  applied 
Dooent  due  to  an  aerodynamic  fbrce. 


ero 


(10a) 


f  -4^(p0t>^  yp) 


CWf  (WeJ 


J 


11. 


(10b) 


‘^aero 


(10c) 


+  CTmGI^D 


““i  Q  AU/y^  -  At^nxpC^ 


Co<.^to^O 


Do  ACPy 

-f-  CmG^D  /CG-CP 


[A(^o6| 


2.  Control  Force  Moment 


+i(C(>M^C»4f  ■t'AU4>^e4i^l9 


This  moment  is  the  result  of  control  forces  obtained  from  either  reaction 


Forces  and  moment  arms  are: 

(11)  Rt  j  5: 


y^Xe  •  ^ 

Axc  " 

Ajffc  •  ^ 

-  -7%-i'h  j 


Thti  control  inO't.ent  Is  defined  as: 


(1^?)  KAc 


BvnluAtlng  each  Indivldviai  croas  product  giveb: 

-  Fx  (-Xc^' -F*h  T' 

-  F.  h  4/  -h  Fi  X. 


^!ibc  K  h  ^  >Cc 

-^jeX  ^  •  fie  Xc  .V 

A  aumatlon  of  the  Individual  control  momante  about  the  body  axes  produces 
the  foilowliiQ  equatlona: 


3. 


d"*)  x.„u  ■  h(-F*+<i+fix-li) 

d*"-)  Sitch*  Xe(F*+(i) 

He„.  .  -Xc(r*  +F*) 

Misaligned  Vector  Mongot 

The  omaent  produced  by  a  oanted,  displaced  thrust  vector  can  be  derived 
as  follows: 


13 


Diagram  2 


Tha  moment  is  expressed  as: 

^  ^ 

where  from  diagram  2;  ^ 

Taking  the  cross  product,  the  moment  becomes: 

(15)  Hth  ■  i 

"“•k.'  |(  Fi-  ■t(F'T 

The  body  axes  cooqponents  of  the  moment  due  to  the  thrust  vector  are  then: 

^roii  -  ^Yr  -  4  Zr  ^ 

**rpitch  - 

^yaw  ■ 


Damping  Moimint 

It  remains  only  to  define  the  aerodynamic  damping  moments  about  the  x'  , 
y'  and  Z*  axes  and  these  are  Mly,  and  Md^*  respectively. 

The  total  applied  moment  is  the  sum  of  the  individual  moments  derived 
previously,  for  equation  6. 

Therefore : 

Body  roll  moment 


(17a) 


iVi  ^  J 

+  ACP,  04-11'  - 

+».  (f>  C4<£j| 

+  h[-r.  +  F.-*-  F»-F^ 
AYr  Ft -  A^r  FV 


Body  pitch  moment 

S  MfiTCh 

(17b) 


-^ACiP* 


•fCai 

-^^Oa^CMisJ 

+  -6.  Hr  FrOt/*^  C**fr-t-  ' 


Body  Yqw  Moment 


7  M 


(17c) 


.  C^QAi^  fC^r-CP\rx  /  J  ^  ^ 

^ 7w(~^///T7  ' 

■^ifj  "*•  <9y 


5..  Total  Rotatloneil  Equationa  of  MotlQ;i 


The  total  rotational  equation  about  the  roll  Jixis  ic  found  by  eqiiatinc 

the  {T.P7  terra)  in  equation  I7a  to  that  in  equation  ';ia  and  dividing  by  roll 

^Pix. 

moiiient  of  inertia. 

^  ^  ^  0 

-P  ^ V«>C  eUa.  04-^-  C^2.0.^-*%.  ^  0  C*X  ^ 

JLlX  ctfZ.  &  Cox.  ^  0 Cffx  ^  —  0Ju>fv  •~^C^a^0.Oon^0.^otx^) 

¥> 

^  F  ^Cgsu0(coxL^  ^F0 

ipbjl 


&  y^-y^  y.^  •*-  ^ 

m 


C^cf>  Cjtrt-  0J 


16. 


r 


Utonli.  <4u>ttng  the  (2/^tch  '•™)  *»  e^uetlon  17b  to  that  in  equation  5h 
and  dlYUtn*  hy  the  pitch  aali  moment  of  inertia,  producen  the  body  pitch 
Axis  rotatlOMl  squat  ion  of  motion. 


~  5ijY^  -‘iiAutoB  - 

(18b) 

1-0^/^  -^I’F^+F^  +  TyUv' 

^j^iC Qo^^idu09M*yF~‘^^^ 

-^tFt  -Kr^rJ^^A-  —  O 

*ir‘  ly' 


rr. 


The  rotational  equation  about  the  body  yaw  axia  Is  derived  by  equating  the 
(3EAf  term)  In  equation  17c  to  that  in  equation  5c  and  dividing  by  the 
yaw  moment  of  inertia. 


(18c) 


-.  Ji 


*  )d*j*Op6^^ 

CbC(h6,^ 

0o>^9  y- 

+•  ^  ■+z^y) /U*ie 

*  +V'^C<Hi*w4Ji/^--d.^en 

'^+>L‘pi+f*7  +7y,^  , 


—  CidGAD/ci^. 

|4u^|  x^\,  l> 


OcMi^  ^ 


/V/  ^  -J 


10. 


I 


S«ctlop  II  -  Translational  Eountlons  o€  Mo r ion 

The  total  translational  equations  of  ci  tionfVrltten  with  respect  to  the 
space  fixed  reference  ayetem  are  now  derived.  The  total  acceleration  along 
any  space  fixed  axle  will  be  due  to  the  num  of  the  aerodynsTilc  force,  the 
weight  component,  the  control  forces,  and  the  thruft  component  alon^  that 
particular  axis,  divided  by  the  maea.  The  force  vector  equations  are: 

r,.  +  f-J 

(19)  Ty  +T.{ 

F,  -  +T*ie 

Aerodynamic  Forces 

By  means  of  a  vector  transformation  the  body  aerodynamic  forces  are 
defined  along  apace  fixed  axes  aa: 


(ao) 


'aero 


P  •  k 


*aero  aero 


Subatltutlng  equations  8  into  equations  20  results  in  the  following: 


(21)  a 


Qfo<.9 
*  Tji^lLiisO  \p 


4i^4i«re^ 


+3 


(21b) 


*  "“0® 


(21c) 


Vfetght 


+^  4^hJ  d  Ai*0'>p(lo^  9 

l^l  J 


(22) 


The  vehicle  weight,  .-jpace  fixed  component!  are:  (Ignoriii,'  trajectory  twiet) 

-4iddJ-Q- 

Th^Y  =  o 


Control  Force 


The  control  force  components  along  'pace  fixed  axes  hre: 


The  thrust  vector  space  fisted  axes  con^neots  are: 

*►  >T  Ce^  y*  •— Gx>^ 

Fry  ■ 

Ft  f-hJdin)^  ^QjdCtO}^ 

Ft  C&^4 

The  Total  Translational  Equations 

The  apace  fixed  accelerations  are  now  presented  as  the  sum  of  components 
of  equat^*^'.  ^1,  ,  23,  and  24  divided  by  the  mass  term. 

X  ^  C<HiOC&^p  — ^^*e»oXl. 


P5») 


-1 


+(ptt^ j Q»-«f  - 

^  ^ 


A  • 

'/'  <5»:y^^Va^^  OfXfPX^^ 

"X  C«a.^^.Xv\^  f£03.^ctx^  XXh.  O^Xts."^) 

X  y4X\/^  (cn.^  .X*}.  ^.Xn  (( ^  ^) 


(a5c) 


2^  **  ^ 


Section  III  -  Approxlnata  Forms  for  Equatlona  of  Motion. 

Preliminary  Slnqpllfl cation 

Although  equations  18  and  23  are  total  equations  and  may  be  programmed 
for  very  large  digital  and  analog  coBq)uter  systems,  they  ere  complex  and 
cunhersGoe  and  can  be  simplified.  The  first  step  In  simplifying  these  equations 
would  be  to  reduce  the  number  of  terms  in  equations  18  by  assuming  that  products 
of  inertia  are  negligible  and  ttiat  for  most  syanetrlcal 

ballistic  missiles.  Further  ir  impli float  ion  is  achieved  by  expressing  body 
axes  moomnt  along  the  Xj  Yj^  Z-j^  axc:j  with  the  following  equations! 


22. 


(26) 


i^Vj  SS 

a  Mv'<W<^  -  M,' 

M 

Incorpormtlng  th*  above  aBBumptlon  and  auhatltutine  the  expreaeiona  for 
»V<  “nd  M,4  Into  equations  26,  yields  the  following  rotational  equations. 

(27«)  ^-V.i^A/aOd 

+1^  ^ 

*  X,'  X«' 


o 


+  It  (^»  ^ 


-JC 
X 


~|c*^4'('i'*'5r)  +  "4«^<^(  fi  +  g  V 

^  Pr  (Ainoff  C®^  ^  .^diiufir 

CW?f  ^W<S^  &K^-  -^^^C®*^d*(8*  ^ 


-X 
J 


VCo^0  '-^yQAi^^ 

+  ^j^CAi>d£iO^  -  -ACPy  fl(Ht^  FrfinJff  Cd4^(U^ 


-I-  CjoC^AD  /co~cP 


xf4U>oclV  I> 


X 


Pcwf  +&)  ->^<<'(fk+&^-4|r  r.Ck«5.(UA>^ 
^  Tr-  (!«4.  Cht^  -  Ft 


+  m<,.G>»4>  ^  '^<*tx^i:,o(b  -  O 


PTturbtioo  Bqwtlon* 


ViM  perturbation  equetions  are  obtained  by  eubtractlng  a  etandard,  eero  angle 
Of  attaek  trejaetory,  froa  equations  (2^)  and  (27) •  Ibe  standard  trajectory  ^an 
be  deeoribed  by  the  following  paraneters: 

Pitch  angle,  ^  ^ ^  ^  steepest 

tr.j.«toru.,  ©;--^ao 

and  there) ore: 

Oi  ir  •  ?e  X  -r  Xg 

04^4  9#  •-4-UuJ^  B  •  ig  X  ■  X, 

2  X  -  Xg 

Control  forces  and  misalignmento  are  zero  aloo. 

Substitution  of  these  psursneters  into  equations  25  and  27  glveu  the  atundard 
trajectory  equations  as: 

x«  tMBt  -h  ql  G^e^-o 

m  cr  TH 

•t 

Y*  «0 

j6*^eg  —  ^  ^ 

V  ••  V 

0.  •  >i  -e,  -o 

Subtracting  equations  (26)  from  equations  25  and  27  give  the  perturbation 
aquations  of  notion  ast  ^assumlnf  OS) 

(2Sk)  «ax  +  c#4.e^) 

■*•  C*iyi)  +  |^(4i«>s 


m 


■A 


+^,(CWv^  -tAi^e4^^ 

+ji 

M 


Ca^y/  ‘fj^4yto(f)At^s  Auwjy^ 

-_  (^j(iw^^  Oh 


•f 

<t>  ~  >-4<:»09 

OtHt-T^  +-4^^  ja^o>^ 

■^1'^/ '^®  -/aH*  Q>*;^ 


^Gyi(p4^y^ 


M 


--^_^.[-F,+  Fa  +  f*''i3 

_  AY^fr  '*’ 


*  Cv  OAO  /4»^]fi.^e  (Wy-t- 
X|AiMj«c|i.  3>  /iJV/  '  n 

_  ^1-  r  (W  4>  (Fa  +  ^)  +  4>  (  Fx  +  fjc^ 

_  P,  Co-«^  + 

+  222^1; <i«H.(li  -  -o 


27- 


(?9f) 


Co^e  ^  AiMj^ 

+  -®S-jAC5;d^4)-ACpy  -^lvC«><|fi^^ 

+  •^rC<H<l.<l>  (»^x+»i:) 

—  ^  Oi^  (S/i  Q^fSy  $iUo^ 

*♦•  — ^  ^  Fj  jU»o^ 

-H  -f-  -32Myi  KuC^^  eaO 


Linearized  Perturbation  Equatlone 


Further  simplification  Is  achieved  vhen  all  angles  except  B  *  XL  and  ^ 


can  be  linearized. 


Incorporating  this  assumption  with  equations  29,  the  linearized 
perturbation  equations  of  tjiot ion,  are  obtained  as; 


(30a) 


+  ^(c»4.e-c<>-4e,9 


(30b) 


-  -^^OKS-Oo^e^- 

—  — ^  ^^^C<hA<^  A*-^e 

-AY  CoA© 

-  -4«4o^  '*■ 

—  .A*/»oe'y'  — .4<^^  *  ^ 


(30  c) 


^if—  *^|>6we->4u*Ufi>^ 

_^-f«^X«)<j>eo^e  -  (-Se^Cik^  CiHie 

+  (AiyQB  ->din}0^  —  £k  /iitOI^  (!,,<^^fi 


— ^  fi/>(!^4>0»<>e  =0 


(30d) 


(30e) 


A«/io©-)^  — 


+  ^C0-4(t(l<HI.©J 
—  C<K,(^ 

+  'j^j'  ^  H"  Ajt^O'Yy 

+  ■‘di^fi  Qo-nt^J 

+  -A  ?r  Fr  tf y  + =0 

r«'  r  ;^  , 


+  '4>^Co^$  Aiioe  + 1** 

- + -■AirtOif)  (  F* -f 

_  A}k./rj^^  ^ 

m3C  mXm  Jm 


;;o. 


t 


(30f) 


y>(U^0  ^  ^  ^ 

+»,  [■^  -  4|fV  <^<H»^  + 

+  4)C^^-*C'4+/ir;;^^ 

-  -^Ft-4^*  +  -^  fyCiMt^  +-^FT?r(W 


P>rlv«tlon  of  th«  coipon«nta  of  angle  of  attack 

Thai  norml  force,  T^,  can  be  vrltten  In  apace  fixed  componenta  by 
expreeainc  along  space  fixed  axes.  Tti»  final  eiqpression 

la  the  following: 


(»)  ^  “  CM,<iA^4  (^^*9  -|^^y-(We  -h |i 

-  ic  +jiy'A;H>edaMr«+^^a,i^ 


Because^  is  aaall,  tbe  coappneni  along  the /(iT  bxI*  i*  approximately  equal  to 
C^^9At^(U»  and  the  J  cooponent  la  approximately  equal  to  C^QA<^y  . 


31. 


mmrn 


I 


Therefore: 


4- *  C^*’6?7 

jV/'^  ®U 

CM,<5<<ioty  -  CNj^GLAj^Ayfl^e-i^i- 


The  Angles  of  Attack  In  pitch  and  yaw  are  defined 
then  as  follows: 


0j^e-(--^y>4^e  +  j^Ca<e 

(32) 

•  •  • 

<^y  -  -j^  -  y'xl^e  Oxe 

The  rates  can  be  expanded  to  a  standard  rate  and  a  perturbation  rate  as: 

X  -  X,  +AX 

Y  »  Ys  4*  AY 

2  »  +  A2 

also 

•j^  »  (WGs  --<dut>^5  - 


Substituting  the  above  Identities  into  equation  32  gives  the-  following: 

(where  (0—0^  and  C^it/A/ )  are  treated  as  snail  quantities 

Just  as  10^) 


{33a) 


m 

a 


M 


32. 


(33b) 


Wind  affacts  can  b«  added  to  the  Qif^end  eqjuetlone.  The  terms  In 
equstlone  30  that  can  be  further  reduced  by  aubatltuting  or  O^V 
where  possible  are: 


Ucom, 

fm=cm.  -  G.«l«|,aA 

'TT) 

(3)  4t4o9  -h  ^  -h  6^1^ 

becomes 


(W© 

•TTl 


33. 


(M 


oecotoes 


(d<:  pCo-a^  +  Ody  4^4)j 


(->) 


*”  4-ji  'h44^4>  ^$f) 


becomes 


(g(pA*^0-  OCy 


C<HL 


(6) 


j]^W>9-4- j^y/^ef 


becomes  Cm^GA 


OTy 


and  Anally, 


(7) 


'^|jw>'-|vr'^ 


j 


becomes 


^^eiGiA 

X 


Substitution  of  the  above  Identities  back  Into  equations  30  presents  the 
linearized  8lx-degree>of-freedo(r  perturbatiun  equations  as; 


(34a) 


(34b) 


((UM4®--4^^^)  •+  Odp>AuiO^ 

—  fF  ~h  =0 

21  y  -f-  ^  Ot»-tf  0  -  0i.y  -  ^ 

*”  — 1^-  do-^s^y/ 


=o 


■I 


O^c) 


-^^j4i*>e-Q**$J  --^fYG<>^^'^4> 
— Qf^e  =o 


0kA) 


«»  ^ 

‘‘^/y^AAO  0 fiJK*B 

-  -fc  [-  ■"  ^•■*  -  ''-3  -  '^ '''''' 


^  «■*  k  J  •  ••"J 

^^Frfiy  +  -^’‘'=0 


(3*^6) 


§  -f.  C<K.«9  Q/>^g 4it0^ 

+  (!»4/SAi^6  +1), 

+  ^^(^)ocr 

-  rr(f^lU^^  ■h^r'^0- 

^  XT 


j  • 


(3'*f) 


y  -h  y^Qjd^^ 

"“CPtt+ 

—  Fr'‘di^^  -h  Fj-Q(K(f)  + 

- ^  Fr^  -/.  <^!i.' y.  /!dic!/l>i40^=0 


i3‘‘t) 


«w>  .  e  4-  F  ^ 

oCy=  yte^^-  j^ 


Section  IV  -  Body  Moment  Derivation  B 

The  body  tnomenta  aa  expressed  In  equations  3  can  be  equated  to  the 
respective  applied  moments  of  equations  17  to  give  the  following  equations: 
Body  Roll  Moment 


(35a) 


“1^ 

4-  ^ 


ACP^ 


-r 

+  h^Fz  +  r^+f^-Fjs^  + 


.>L>. 


Body  Pitch  Mo—nt 


(35b) 


Body  Yav 
(35c) 


*  -A  iicg  (^fi  +  /i]  +  ^*T  Ft  Co*  Pp 

« 

--4uw/  fi^tf 


VbmvX 

»  Dl,ACI^-Xe  [r^-^rJ  -Xrr^OK^jUtfiy 

9 

+  -J^  >4mi>  ^ /C^  ® 

—  AYtFt  C»<#Y  —  M<i*' 


» 


37. 


When  the  corresponding  standard  trajectory  equations  are  subtracted  from 


those  in  eqxiatlons  35>  the  perturbation  equations  of  the  body  moments  are 
obtained. 

The  perturbation  equations  of  body  moments  are  simplified  by  incorporating 
the  aame  assumptions  as  in  derivation  A. 

The  resulting  perturbation  equations  are: 

Body  roll  moment 

-  (ocp 

+  h[-Fx+F3,+li.-1^] 

~  ^^rFr/fy  —  Mdy! 


(36a) 


-fix  =  - 


Body  pitch  moment 

(36b)  *  ~X>,AC^  +Xc[F„-(-y 

—  Cutset Al>  +  O^y 

+•  FrZiZr  +-  x-tFt^p  -MAy' 


Yaw  Moment 


(36c  ^ 


ri,-  +-^^(x/-xj=  _x,[Fi+r^j 

— NAcAg^ 

"i”  '^j^OCp  ^/Li>»j(^  (j^ 


wtwre  ^f,  and  have  been  defined  pn;viou8ly  .xis: 


^A^Ckhde 


yj  Co-a^  — 


|Vl=  ^Vi7p 


w. 


References 


1.  Airplane  Performance  Stability  and  Control  by  Courtland  D«  Perkins 

and  Robert  E.  Hage  published  by  John  Wiley  &  Sons,  Inc.  copyright  19**'9. 

2.  Ballistics  of  the  Future  by  Dr.  J.M.J.  Kooy  and  Prof.  IR  J.U.H. 
Uytenbogaart  published  by  McGraw  Hill  Book  Coriipany,  Inc.  copyright  19^. 

3.  AME-TN-7-60  Non-Linear  Six-Degree  of  Freedom  Equations  of  Motion  for 
a  Ballistic  Missile  by  R.J.  Thompson  and  G.A.  Socks,  l6  March  i960. 


Nome DC la cure 


Symbol 

A 


rc 

T,  Z) 

« 

h 


V 


I*' 

J3&' 

Jxy' 

JyV 


Definition 

Missile  reference  area 

Gradient  of  the  nornial  force  coefficient 

Normal  force  coefficient 

Axial  drag  term 

Reference  diameter  of  missile 

Aerodynamic  normal  force 

Total  aerodynamic  force 

Thrust  of  each  control  meiiber 

Engine  thrust  vector 

Force  with  coiig>03ents  in  X,  Y,  Z  space  axes 
system 

Axial  dislence  between  the  center  of  gravity 
and  the  center  of  pressure  of  the  body 

Acceleration  due  to  gravity 

Distance  of  control  force  from  body  fixed 
longitudinal  axis 

Moment  of  niouentum  about  body  roll  axis 
Moment  of  momentum  about  body  pitch  axis 
Moment  of  momentum  about  body  yaw  axis 
Moment  of  inertia  about  body  roll  axis 
Moment  of  inertia  about  body  pitch  axis 
Moment  of  inertia  about  body  yaw  axis 
Product  of  inertia  about  body  pitch  axle 
Product  of  inertia  about  body  yaw  axis 
Product  of  inertia  about  body  roll  axis 
Body  roll  rate 


41. 


Units 

1/RAD 

4g 

M 

Kg 

Kg 

Kg 

Kg 

Kg 

M 

M/sec^ 

M 

Kg  M  sec 
Kg  M  sec 
Kg  M  sec 
Kg  M  aec^ 

w  2 

Kg  M  sec 
Kg  M  sec^ 
Kg  M  sec^ 
Kg  M.sec^ 
Kg  M  sec^ 
rad/sec 


i 


q 

r 


M 


^ero 


Mrh 


to 


Mix 

Mdy 

Mdz 

4 


Ro 

R 

Rin 


r,j, 


I 

J 

k 

k' 

SI 


Body  pitch  rate  rad/eec 

Body  yaw  rate  rad/eec 

Moment  Kg  M 

Total  aerodynamic  body  moment  Kg  M 

Moment  due  to  the  control  aystcm  Kg  M 

Moment  due  to  thrust  vector  misalignment  Kg  M 

MlsBlle  weight  vector  Kg 

Mass  of  the  missile  Kg  sec^/M 

Damping  :i<ODient  along  body  roll  axis  Kg  M 

Damping  moment  along  body  pitch  axis  Kg  M 

Dancing  moment  along  body  yaw  axis  Kg  M 

2 

Dynamic  pressure  Kg/M 

Mean  radius  of  earth  M 

Distance  from  the  center  of  the  earth  to  the 
missile  CG.  M 

Initial  distance  from  the  center  of  the  earth  M 

to  the  missile  CO. 

Distance  from  body  center  of  gravity  to  ap¬ 
plication  of  thrust  vector  M 


Unit  vector  along  space  fixed  X  axis 

Unit  vector  along  space  fixed  Y  axis 

Unit  vector  along  space  fixed  Z  axis 

Unit  vector  along  body  fixed  X  axis 

Unit  vector  along  body  fixed  Y  axis 

Unit  vectu-  along  body  fixed  Z  axis 

EJulerian  ntigles  of  rotation  of  the  body  axes 
around  the  X,  Y,  Z  space  fixed  axes  system  rad 

Angle  at  the  earth's  center  between  initial 

local  vertical  Rnd  the  local  vertical  at  any  time  rod 
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Thrust  vector  mlMllgQir«Dt  ancles  (see  diagram  2)  rad. 


Standard  trajectory  Euler  pitch  angle  rad. 

Total  angle  of  ett«c)<  rad. 

Co-Tiponents  of  oC  in  the  pitch  and  yaw  planes 
respectively  rad. 

Angle  between  the  original  local  vertical  and 

the  velocity  vector  rad. 

Distance  from  body  ('.C.  to  application  of  thrust 
vector  parallel  to  the  longitudinal  body  axis.  M 

Distance  along  the  body  longitudinal  axis  to 

the  control  force  station  from  the  body  CO.  M 

Missile  velocity  H/aec 

Space  fixed  standard  trajectory  velocities  M/sec 


Body  Eulerlan  axes  system 

Space  fixed  axes  system 

Dlsplacenent  of  the  thrust  vector  In  the 

body  Y^and  z' directions  respectively  (see 

diagram  2)  M 

Displacement  of  the  center  of  pressure  on 

the  Y  axis  M 
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Displacement  of  the  center  of  pressure  on 

the  Z  axis  M 
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Momenti  of  Momeatum  Referred  to  Moving  Axi» 
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